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The description is given for the speckle interferometer of the BTA 6-m telescope of the SAO 
RAS based on a new detector with an electron multiplication CCD. The main components 
of the instrument are microscope objectives, interference filters and atmospheric dispersion 
correction prisms. The PhotonMAX-512B CCD camera using a back-illuminated CCD97 
allows up to 20 speckle images (with 512x512 pix resolution) per second storage on the 
hard drive. Due to high quantum efficiency (93% in the maximum at 550 nm), and high 
transmission of its optical elements, the new camera can be used for diffraction- limited 
(0.02") image reconstruction of 15™ stars under good seeing conditions. The main advantages 
of the new system over the previous generation BTA speckle interferometer are examined. 



1. INTRODUCTION 

Improvement of the angular resolution of 
ground-based optical telescopes remains one of 
the most important astronomical problems. The 
advantages of the largest astronomical instru- 
ments cannot be implemented without an intro- 
duction of the newest methods of correction of 
phase degradations, arising during the light prop- 
agation through the turbulent atmosphere. This 
problem obtains a special meaning in connection 
with the development of the new generation of 
giant telescopes. 

The adaptive optics and speckle interferom- 
etry allow approaching the diffraction limit of 
the telescope's angular resolution. Image recon- 



struction in speckle interferometry is achieved by 
the integration of a series of short-exposure im- 
ages of the object (exposure time is 5-40 ms) 
with the subsequent calculation of the ensemble- 
averaged power spectrum of the object [1], and 
its phase defined from the bispectrum [2]. In 
contrast to the adaptive optics, which is effective 
today mainly in the infrared, speckle interferom- 
etry can be used for observations in visible and 
near UV bands. In addition, speckle interferom- 
etry is realizable under poor atmospheric con- 
ditions, while the adaptive optics always needs 
the best seeing. Another important fact is the 
value of speckle interferometric equipment, which 
is significantly lower than the expenses for the 
development of adaptive optics for a large tele- 
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scope. 

The speckle interferometry technique has 
been significantly upgraded throughout its 35 
years long history. First detectors based on im- 
age intensifiers and film cameras [3] were replaced 
by the television photon counting systems [4]. In 
the 1990s, television tubes were driven out by 
fast CCDs [5, 6], however, different types of im- 
age intensifiers were still used as brightness am- 
plifiers. Lately, these systems were replaced by 
new EMCCDs with electron charge multiplica- 
tion. They have certain advantages for applica- 
tion in the optical speckle interferometry: pho- 
ton counting sensitivity allied with the maximum 
quantum efficiency, fast image readout, and high 
geometric and photometric stability. 

In the present paper we give the description of 
the EMCCD-based speckle interferometer (SI), 
which is in active use in the observations con- 
ducted at the BTA 6-m telescope of the Spe- 
cial Astrophysical Observatory of the Russian 
Academy of Sciences (the SAO RAS) since 2007. 
The main characteristics of the optical compo- 
nents of the instrument are given in Chapter 2. 
Then, Chapter 3 describes the EMCCD detec- 
tor. The control system is briefly presented in 
Chapter 4. In conclusion, capabilities of the new 
instrument, estimated from the observations of 
binary stars, are laid out in Chapter 5. The al- 
gorithms and software for image reconstruction 
from the series of short-exposure speckle inter- 
ferograms will be described in a separate paper. 



2. MAIN OPTICAL COMPONENTS 

The main optical elements of the SI optical- 
mechanical unit are (Fig. 1): the electro- 
mechanical shutter with the 6 mm aperture to 
block the light beam from the telescope (1), mi- 
croscope objectives for the image magnification 
(3), the set of prisms for atmospheric dispersion 
compensation (4), a the set of interference fil- 
ters for the narrow spectral band selection (5). 
The /:4 beam from the telescope's primary mir- 
ror forms an image in the focal plane (2), and 
then it is recorded by the EMCCD detector af- 
ter the passage through all optical elements (6). 
Let us consider the main components of the SI 
in details. 

2.1. Microscope Objectives 

The SI microscope objectives are designed for 
matching the image scale in the prime focus of 
the BTA with the detector's pixel size. The size 
of a single speckle in a speckle interferogram, cor- 
responding to the diameter of the first dark ring 
in the Airy disk at the wavelength A, is equal to 
d = 1.22A/D, where D is the aperture diameter. 
In the BTA focal plane, for Ao=550 nm d=2 fj,m. 
The detector pixel size is 16 fj,m, therefore, fol- 
lowing the Kotelnikov's theorem {fdiscr > "^fmax, 
where fmax is the maximum frequency in the sig- 
nal spectra), an image magnification of approxi- 
mately 20 times must be applied to reconstruct 
the speckle profile. 



Figure 1. Components of the mechanical optics unit of tlic BTA speckle interferometer: 1 — shutter, 

3— changeable microscope objectives, 4— atmospheric chromatism compensation prisms, 5— a set of 
interference filters , 6 — EMCCD. Focal plane of the telescope is enumerated 2. Black arrows point in 

direction of the components movement. 
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Figure 2. Spectral characteristics of interference filters. The transmission in maximum and half widths (in 

brackets) are marked over each curve. 



We use two high-quality Carl Zeiss objectives a magnified, inversed and free from aberrations 

with the corrected field curvature to match the image of the star at the CCD. With the BTA 

scales: a 2.5-fold achromat with the numerical main mirror curvature radius of i?=48051ibl5 

aperture 0.08 and a 16-fold plan-achromat with mm (floor measurement by L.I. Snezhko using 

the numerical aperture 0.35. They both have the Hartmann method), the microscope objec- 

a back working distance of 160 mm and form fives provide scales of 3.44"/mm and 0.54"/mm 
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on the detector. For the 512x512 pixel detec- 
tor with the size of the sensitive area 8.19x8.19 
mm, the corresponding scales are 0.0550"/pix 
and 0.0087'7pix, giving the fields of view of 28.2" 
and 4.4". The 16-fold objective is the main work- 
ing objective of the SI, while the small magnifi- 
cation is used only for pointing on the object and 
tying-in the position angles of the system. 

For some purposes that require obtaining a 
detailed profile of the Airy disk, one might need 
to use higher optical magnifications. One of the 
examples is the image reconstruction of visible 
disks of cool supergiant stars. For such cases, 
we anticipated an installation of a 32-fold micro- 
scope objective with a corresponding reduction 
of the field of view. 

Small convergence of the beam after the 16- 
fold objective (1:64) permits fitting other optical 
elements (prisms, filters) on the axis between the 
objective and the detector, without resorting to 
additional transfer optics. 

2.2. Interference Filters 

To keep the coherence during the interfero- 
metric picture recording, the width of the filter 
transmission band AA is selected following the 
relation 



where A is the center wavelength of the filter 
transmission band, ro is the atmospheric coher- 
ence radius (Fried parameter) [7], or 



A (Jr 

where cr^ is the mean-square aberration of the 
incoming wavefront phase given in radians. 

Under the seeing of approximately 1", (ro~20 
cm), the filter limitations for the BTA applica- 
tions are chosen by the ratio A A/ A Ril/30. For 
poor atmospheric conditions the filter passband 
must be reduced. However, in most of the cases 
the filter transmission band can be doubled with- 
out significant dissipation of the high-resolution 
data. This simplification leads to a larger de- 
tected flux. In the near infrared, the filters with 
a band width of 100 nm can be used. Narrower 
filters can be installed in the beam to cut out the 
selected bands in the spectrum of an object. One 
of the examples of such tasks is the Mira stars im- 
age reconstruction in the TiO absorption bands 
and in the nearby continuum. 

Five interference filters, manufactured by An- 
dover (USA), are arranged on the turret assem- 
bly, and are being brought in the beam by the 
rotating mechanism. All the filters have the di- 
ameter of 25 mm (the light diameter is 21.2 mm) 
and thickness of 6 mm. The parameters of the fil- 
ters are given in Fig. 2. The thickness of the filter 
dielectric layers is constant (0.25%) throughout 
the surface in order to ensure the design perfor- 
mance. During the detector calibration proce- 
dures or the system testing in the white light, 
no filter is needed in the optical path. For these 
purposes the filter wheel has one additional hole. 
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Figure 3. Dependance of the atmospheric 
dispersion from the zenith distance for three 
wavelengths. 

2.3. Atmospheric Dispersion Correction 

The atmospheric dispersion has an effect on 
the resolution in speckle interferometry even 
while applying narrow interference bands. The 
dispersion value can be estimated from the em- 
pirical formula proposed by Lambert as far back 
as in 1759: 



dz 
dX 



tan z dn„ 



dX ' 



(3) 



where Ua is the refraction index of the atmo- 
sphere, z is the zenith distance of the object. 

The Owens empirical expression [8] is most 
frequently used to describe the wavelength de- 
pendence of the refraction index 
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(130 - A-2) 
6487.31 + 58.058A"2 (4) 



(38.9- A-2)_ 

-0.71150A-^ + 0.08851A-^ 



where Dg and Dy^ are the factors accounting for 
the density of the dry air and water vapors, re- 
spectively. The value of ^^(A) can be estimated 
with an accuracy of 10^^-10^' from formula (4). 
Fig. 3 shows the dependence of zenith distance 
from the refraction index in three wavelengths 
for the BTA site (pressure 720 mbar, tempera- 
ture 0°C, relative humidity 30%) as follows from 
expressions (3) and (4). 

For example, an image of a point source at 
z=30° and A=500 nm will be shifted by 0.05" 
from its true position. The vertical extension 
of speckles in the region A=500 nm through the 
AA=20 nm filter will be around 0.02", which is 
comparable with the size of the diffraction spot. 

The simplest way to compensate the atmo- 
spheric chromatism is the use of a direct-vision 
prism with the dispersion vector opposite to the 
atmospheric dispersion. The angular dispersion 
of the prism must change with the zenith angle. 
One should keep in mind that in the SI we use the 
16-fold magnification, therefore the prism correc- 
tor, which is installed behind the microscope ob- 
jective, must provide a proportionally higher dis- 
persion. 

Beam deviation a in the prism with the index 
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of refraction n and the apex angle (3 is equal to 



a = (n(A) - 



(5) 



For a prism consisting of a pair of wedges with 
the angles /3i and (32 and the indices of refraction 
ni and n2, the refracting angle is 



(ni(A) - - (n2(A) - l)p2. 



(6) 



The wedge angles ratio is obtained by speci- 
fying the zero deviation at Ao=550 nm: 



ni(Ao 



1 



(3i n2(Ao) - 1 



(7) 



To compensate the atmospheric dispersion at 
the maximum zenith angle 2=60° in the 20 nm 
band at Aq =550 nm, the following equation 
should be valid: 



dz 
dX 



dni 
~d\ 



132 



dn2 
~d\ 



(8) 



leading to the apex angle of the first wedge: 



(9) 



-n2(Ao)-l 

where A means the variation range of z and n in 
the wavelength range from Aq — 10 nm to Ao + 10 
nm. To correct the atmospheric chromatism at 
smaller z, the compensating prism should pro- 
vide lower dispersion. 

In our observations until September 2007, a 
direct- vision Amici prism from a pair of wedges 
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Figure 4. Dispersion curves of the prism (solid 
line) and the atmosphere (dashed line). 



has been used for dispersion compensation. The 
dispersion value has been changed by transla- 
tion of the prism along the axis: the lowest value 
was obtained in the immediate proximity of the 
detector; the highest dispersion was reached at 
the remote position. The advantage of such con- 
struction consists in its low Ught losses in the 
glass and on the surfaces. The disadvantage of 
the assembly is its remaining dispersion in the 
circumzenithal zone. This dispersion cannot be 
eliminated as the prism cannot be installed closer 
than some minimum distance from the image ac- 
quisition plane. 

The dispersion curve for the prism compared 
to the atmospheric dispersion is shown in Fig. 4. 
It can be seen that by the selection of the wedge 
angles and refractive indices, it is possible to 
match the dispersion variation of the prism and 
atmosphere in the range of 550-700 nm. 
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The new compensator was built on the base of 
a direct- vision Risley prism consisting of a pair of 
identical prisms rotating around the instrument 
axis in opposite directions (see Fig. 1). The com- 
pensator gives a dispersion of 27sin(0/2), where 
7 is the dispersion of a single prism, 6 is the ro- 
tation angle of the prisms. Each prism in the as- 
sembly is a glued pair of wedges (crown -|- flint) 
with carefully selected indices of refraction, pro- 
viding a dispersion of 330"/20 nm at the direct- 
vision wavelength of 550 nm. Shorter and longer 
wavelengths deviate symmetrically in opposite 
directions from the 550 nm beam. The net de- 
viation of the beam is the vector sum of the two 
separate deviation vectors. If the prisms are ro- 
tated by 90°, they combine to act as a parallel 
plate with no net angular deviation. To simplify 
the construction, the rotation mechanism of the 
prisms is made dependent, providing the rotation 
of the prisms by equal angles about the optical 
axis. The rotation angle of the prisms is defined 
by the computer when the object is chosen and 
its zenith distance is calculated. The calibration 
of the 0{z) relationship is made using the obser- 
vations of bright stars at different z during the 
nights of good seeing. 

Due to different distances of the prisms from 
the detector plane (the front prism is at larger 
distance and produces larger dispersion), the net 
deviation of the assembly is shifted from the ver- 
tical. The value of this shift is zenith angle 
dependent. The remaining residual dispersion 
can be easily compensated via rotating the prism 



mounting by a small angle (about 10°) in the op- 
posite to the differential dispersion direction. 

The Risley compensator, consisting of a pair 
of rotatable wedged elements, provides the cor- 
rection of the atmospheric chromatism of speckle 
interferograms in the range of zenith distances 
from 0° to 60°. The prisms were manufactured 
by B.Halle Nachfolger GmbH, Germany, from the 
SF L6 crown glass and heavy flint glass LaSF 
N30. The apex angles of the wedges made of 
crown glass and heavy flint glass are /?i=41.6° 
and /?2=41.9°, respectively. The prisms have a 5- 
layer antireflection coating for the range of 450- 
1100 nm. 

3. DETECTOR 

3.1. General Characteristics 

The PhotonMAX-512B camera developed by 
Princeton Instruments, USA, on the base of elec- 
tron multiplication CCD97 (EMCCD) is used 
as the detector in the SI. The back-illuminated 
CCD has 16x16 fim pixels in a 512x512 frame 
transfer format with the total photosensitive area 
of 8.19x8.19 mm. The quantum efficiency of the 
CCD is QE=93% at 550 nm. The quantum effi- 
ciency curve at different wavelengths is shown in 
Fig. 5. 

The EMCCD has two independent readout 
registers on the output. The first one is used 
for standard, high dynamic range application as 
in traditional CCDs, the second one is devoted 
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Figure 5. Rated curve of the CCD97 quantum efficiency. Dashed line marks the wavelength corresponding 

to the maximum quantum efficiency of 93%. 



for high-speed, low hght level multiplication gain 
mode. In the normal readout application (no 
multiplication), the readout mean-square noise is 
45 e~ and 60 e~ at operating readout frequencies 
5 and 10 MHz, respectively. Under maximum 
multiplication of the charge generated in the pix- 
els, the readout noise falls below 1 e~. The Pho- 
tonMAX can offer a multiplication gain factor 
between xl and xlOOO. The typical relationship 
between the 12-bit setting of the DAC (digital- 
to-analog converter), which defines the voltage 
applied to multiplication register clocks, and the 
charge multiplication factor is close to exponen- 
tial (see Fig. 6). The maximum multiplication is 
achieved under the voltage of about 43 V on the 
extended multiplication register. Further regis- 
ter voltage gain can destroy the CCD. Most of 
the time we use the PhotonMAX-512B camera 
in the high-multiplication mode, and the normal 
mode (no gain) with low multiplication is used 
only during the speckle observations of bright 



stars. 

To reduce the dark current, the EMC CD is 
cooled by a 4-stage thermoelectric Peltier cooler. 
It is stabilized to within ±0.05 °C by a feedback 
control system. The lowest achievable temper- 
ature of the CCD is -70°C providing the dark 
current below 0.01 e~/s per pixel. 

The CCD is mounted on the cold finger in the 
vacuum module. The photosensitive area of the 
CCD is protected by a single optical entrance 
window, which provides low light losses. The 
heat produced by the Peltier device is removed 
from the camera by an inbuilt fan. The fan is de- 
signed for low vibration and does not adversely 
affect the image. 

The camera can be used in the frame trans- 
fer mode in the wide range of exposures between 
10 jUS to 10 min. Maximum readout rate of the 
camera in the full-frame mode is 29 frames per 
second (fps). As the PhotonMAX camera uses a 
frame-transfer CCD, it does not incorporate any 
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Figure 6. CCD97 charge multiplication factor 
dependence on the code setting at the input of 
digital-to-analog converter. 

electro-mechanical shutters that set the exposure 
time. The exposure time is set electronically by 
the software. The external shutter 1 (see Fig. 1) 
is used solely to block the incoming light in the 
SI. 



fer function of the atmosphere is not stable for 
longer periods. To check the time stability of 
the detector, we have accumulated a few sets 
of 1000 uniform (flat) light fields obtained at 10 
MHz operating frequency with an exposure time 
of 10 ms. In Fig. 7 presents the variations of tlie 
summed signal from the CCD. It can be seen that 
the CCD signal intensity varies exponentially by 
almost 10% during the acquisition of 1000 frames 
(54 s). This circumstance must be taken into ac- 
count for photometric measurements. To stabi- 
lize the multiplication level of the CCD before 
the main exposure, the CCD has to be read for 
more than 5 min. 



3.3. Readout Speed 



3.2. Time Stability 

Basic properties of the PhotonMAX camera 
were studied at the laboratory and on the tele- 
scope. The study included measurements of the 
wavelength dependence of the quantum efficiency 
of the SI, an analysis of the system's stability, its 
frame rate in different acquisition modes, the sig- 
nal homogeneity over the field, and the flat field 
response of the system. 

Speckle interferometry operates with series of 
a few thousand short-exposure images of an ob- 
ject. The total integration time of the series 
is less than a few minutes as the optical trans- 



Regardless of the exposure time, the highest 
readout speed of the 512x512 pix full format im- 
age is 29 fps. Binning (combining pixels into 
one super pixel) and sub-region selection allow 
increasing the sensitivity and frame rate of the 
PhotonMAX-512B. The camera allows the bin- 
ning of 2x2, 4x4 and 6x6 along either direc- 
tion of the CCD. The maximum frame rate vs. 
binning for different readout regions is shown in 
Fig. 8. The maximum rate of 350 fps is obtained 
for 6x6 binning and reading only the 1/64 area 
of the CCD. However, all types of binning are 
rarely used in speckle interferometry as it reduces 
spatial resolution. 
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Figure 7. CCD flat fleld time stability. 
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Figure 8. Frame rate dependence of the CCD reading area without binning, and witli 2x2, 4x4 and 6x6 
binning. The dimensions of the reading zone area are protracted on the X-axis. 



3.4. Flat Field Signal and the Dark-Charge 
Pattern 

Even under uniform incoming illumination 
(flat field) of the CCD, the obtained image is not 
perfectly flat. Inhomogeneity of the response is 
explained by different reasons: 

- lack of uniformity of the chip; 

- contamination and granularity of the cam- 



era's entrance window; 

- imperfections of the CCD readout elec- 
tronics, which become apparent mostly on 
the edge zones of the frame. 

Flat field problems do not play such an impor- 
tant role in speckle image reduction as compared 
to the classical long-exposure astronomical imag- 
ing due to the fact that interferometry exploits 
ensemble-averaged power spectra. Nevertheless, 
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the standard reduction procedure in speckle in- 
terferometry also implies a flat-field reduction of 
each speckle interferogram. 

We have studied the SI flat fleld by integra- 
tion of the twilight sky. One of the examples 
obtained by averaging of 2000 short exposures is 
presented in Fig. 9. The fleld inhomogeneity (sig- 
nal dispersion) for this image is 0.5%. Summing a 
larger number of images allows determining the 
PhotonMAX-512B flat flled with the noise dis- 
persion under 0.05%. 

The CCD ofl'set signal measurements with no 
light coming into the camera (bias) were per- 
formed with the closed shutter using 10 /v,s ex- 
posures. An average of 500 frames obtained 
with the highest CCD multiplication is shown 
in Fig. 10. The general view of the base plate 
is shown on the left, and its 3-dimensional pre- 
sentation is on the right. Pixel values of the 
charge are given in the corners of the bit map. 
It can be seen that flrst lines in the frame have 
significantly higher charges due to the charge in- 
leakage from the chip edges. Maximum differ- 
ence in the counts is 0.6% from the full range 
of the analog-to-digital converter scale (65536). 
The bias pattern is fully subtractable if the dark- 
charge pattern is acquired under conditions iden- 
tical to those used to collect the actual speckle 
images. 



3.5. Single Photon Sensitivity 

During the speckle observations of faint stars 
with the maximum gain, the PhotonMAX-512B 
acquires images in the photon counting mode. 
Fig. 11 shows the 20 ms exposure image of 
LkHA 198 (mv'=14.3™) obtained at the 6-m 
telescope in the 600/40 nm filter under 1" see- 
ing. The cross-section illustrates intensity max- 
ima corresponding to single photon events. The 
CCD noise dispersion for this image is 30-40 
counts. Individual photons are detected by the 
PhotonMAX-512B with a signal-to-noise ratio 
of 10 to 50. 

4. CONTROL SYSTEM 

The SI control system diagram is given in 
Fig. 12. Control of optical elements is performed 
by the stepper motor drive controller TMCM- 
310 from Trinamic Motion Control GmbH & Co., 
Germany. It can drive up to three bipolar step- 
per motors through the TMC246 drivers with a 
peak coil current of up to 1.5 A. A built-in micro- 
processor (fJ-P) and a nonvolatile 16 KB memory 
(EEPROM) with the aid of the parameter set- 
ting module TMC428 allow fiexible adjustment 
of the drive for selected translation mechanisms. 
The controller is plugged in the COM-port of 
the control computer through the standard se- 
ries RS-232 interface via the 8-channel fiber optic 
LWL converters. 

The microscope objectives are translated by 
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Figure 9. Image of a small area of the twilight sky integrated in sets of 2000 frames (a fiat field). A cross 
section along one of the lines in the bottom part of the frame is shown on the right. 




Figure 10. View of the CCD background and its 3-D presentation. 



the precision linear motorized translation stage the beam with the positional accuracy of 0.01°. 
ERLIC 85 produced by OWIS GmbH, Germany. The atmospheric dispersion compensation prisms 
The setting accuracy is 2 /xm for the 50 mm are rotated by stepper motors through a differ- 
travel. ential mechanism. Special guidance provides a 

practically slip-stick free movement as well as 



The filter wheel is fixed on the rotation stage 
mechanism DMT-65. A filter can be installed in 



high load capacity for all translation stages. 
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Figure 11. Left: a speckle image of the star LkHA 198, acquired by the PhotonMAX-512B in the maximal 
CCD signal multiplication mode. Top right: a magnified segment with the maximal flux level. Bottom right: 
an arbitrary cross-section, demonstrating the amplitudes of independent photon events. 



The camera and other units of the SI are con- 
trolled by one computer based on Intel Pentium 
4 3.2GHz processor. We use the Shuttle XPC 
SB81P Barebone system with low electromag- 
netic radiation. The PhotonMAX-512B camera 
control is conducted by the ST- 133 controller 
through a PCI-bus via the fiber optic AS-0230 
receiver-transmitters. RAM is increased to 2 
GB to enable rapid accumulation of up to 1940 
speckle images, which are then recorded on the 
computer hard drives in the binary files of SPE 
type with a 4100 byte header. The file header 
contains all the important camera settings that 
were used during the accumulation of the speckle 
images. The accumulation of the data can be 
performed on the high-capacity external hard 



drives via a USB-interface. 

During the observations we often need to 
monitor the accumulation of power spectrum of 
the object in the real time mode. The infor- 
mation on duplicity (multiplicity) of the object 
at the initial stage allows assessing the appro- 
priateness of further acquisition. To solve this 
problem we are currently implementing recon- 
figurable multiprocessor computing systems into 
the process of observations and data reduction. 
Unlike systems based on rigid architecture, re- 
configurable systems allow on-tlie-fly changes of 
the calculation parameters. As a result, users 
may fit the system's architecture to the structure 
of a given task. 

To calculate the power spectrum in real time 
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Figure 12. Block diagram of the SI control system. 
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Figure 13. SI quantum efficiency corrected for light losses on the optical elements: 1— CCD97, 2— CCD97 

+ PhotonMAX-512B entrance window, 3 — CCD97 + PhotonMAX-512B entrance window + dispersion 
compensation prism, 4— CCD97 + PhotonMAX-512B entrance window + dispersion compensation prism + 
filters, 5 — total efficiency corrected for losses on the miroscope objectives. 
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Figure 14. Individual speckle images of DF Tau binary star. Left: an image, acquired with the CCD 
image intensifier system. Right: image obtained with the EMCCD. 




Figure 15. DF Tau power spectra. Left: the spectrum accumulated on 2000 speckle interferograms with 
the CCD + image intensifier system. Right: the same accumulation with the EMCCD. The arrow points at 

the direction of zenith. 



mode, a unit based on a reconfigurable acceler- 
ator of personal computers (RAPC-25) with 25 
gigaflops processing power [9] was designed in the 
Research Institute of Multiprocessor Computing 
Systems, Taganrog, Russia. This unit imple- 
ments parallel algorithms to correct the optical 
images for atmospheric distortion, allows real- 
time reduction of speckle interferograms, and 
on-the-fly calculations of the power spectrum of 



the observed object's image. The unit is repro- 
grammable, an option that can be used, e.g., for 
object phase reconstruction. 

5. EFFICIENCY OF THE SYSTEM 

Quantum efficiency of the SI was mea- 
sured in 9 wavelengths using the MDR- 
41 grating monochromator, and the OL 
730D programmable DSP-radiometer, 
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which is sensitive in the range of 
200-1000 nm. Fig. 13 shows the system's 
efficiency curves vs. wavelength, adjusted for 
losses at the optical components: the upper 
curve (1) gives the QE of the CCD97 chip, the 
lowest curve (5) presents the total QE of the SI. 
Note that the CCD itself has a 20% efficiency at 
A=950 nm. 

We compared the PhotonMAX-512B camera 
performance with the previous detector, which 
has been used for speckle observations until 2007. 
This camera was built on the base of a 3-stage 
electrostatic-focusing image intensifier AEG1510 
(Siemens, Germany) optically coupled to the fast 
CCD camera SensiCam 370LL, manufactured by 
the PCO Computer Optics, Germany. 

Individual speckle interferograms of the bi- 
nary star DF Tau {my=12"^) obtained with 
the 6-m telescope in October 2001 (intensifier- 
SensiCam) and November 2007 (PhotonMAX- 
512B) are shown in Fig. 14. Both observa- 
tions were performed under l"-1.5" seeing in the 
A/AA=800/100 nm filter. In 2001, a scale of 
0.004"/pix was provided by the 32x microscope 
objective, while in 2007 a scale of 0.0087"/pix 
was provided by the 16 x image magnification. 
The difference between the two images is evident. 
Photon noise is dominating at the intensifier- 
SensiCam speckle image. To reconstruct an im- 
age of the binary from such data, it is necessary 
to process up to 10^ individual speckle interfero- 
grams. On the other side, the speckle structure 
is clearly seen on the PhotonMAX-512B image; 



in this case a diffraction-limited reconstruction of 
the DF Tau image can be obtained from several 
tens of frames. 

In speckle interferometry, the power spectrum 
of the star's image is recovered from a series of 
short-exposure speckle interferograms using the 
Labeyrie's method [1]. For binary and multi- 
ple star speckle interferometry at the BTA tele- 
scope we use the procedures described in [6, 10]. 
Fig. 15 represents power spectra for DF Tau av- 
eraged from 2000 speckle images, recorded by 
the intensifier-SensiCam (left) and PliotonMAX- 
512B (right) cameras. The diffraction cutoff fre- 
quency of the aperture in the used filter is marked 
by circles. The period of fringes in the power 
spectrum corresponds to the angular separation 
between the components of the binary (about 
0.1"). Orientation of the fringes determines the 
relative position angle of the binary components; 
the magnitude difference between the two stars 
can be found from the fringe contrast. 

The drop of intensity with increasing spatial 
frequency on the left image is caused by an ad- 
ditional photon noise spectrum, resulting in a 
frequency-dependent bias in the fringe contrast. 
Correcting this bias poses a separate intricate 
problem, being an obstacle for differential pho- 
tometry of the components of faint sources. In 
the power spectra accumulated on the EMCCD 
data, the influence of photon noise leads to a 
common change in the level of background, which 
is easy to deflne. High contrast fringes are indica- 
tive of the power spectrum of DF Tau, obtained 
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in 2007 with the new system. These fringes can 
be traced up to the cutoff frequency, determined 
by the diffraction at the telescope aperture. 

6. CONCLUSIONS 

In conckision, we shall list the main advan- 
tages of our new system based on the EMCCD, 
as compared to the previous generation BTA 
speckle interferometer based on an image inten- 
sifier optically coupled to a CCD camera: 

- reduction of the total object exposure time 
thanks to the 4-fold performance advan- 
tage of the EMCCD; 

- improvement of limiting magnitude up to 
15"* owing to the 10-fold quantum efh- 
ciency advantage of EMCCD over the im- 
age intensifier photocathode; 

- amelioration of image reconstruction accu- 
racy thanks to the absence of photon bias, 
which is the main problem of image inten- 
sifiers application; 

- increase of the maximal achievable magni- 
tude difference in binary star observations 
up to 5"* owing to the high dynamic range, 
and high signal-to-noise ratio in the recon- 
structed image; 

- geometrical stability of the EMCCD-based 



detector and lack of distortions inherent to 
image intensifiers; 
- broadening of the speckle interferome- 
ter's working field up to 4" owing to the 
smaller optical magnification provided by 
the small EMCCD pixel size (16 fj,m). 
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